LU TP 97/03 
NOR, DITA-97/22 N/P 
hep-ph/9704212| 
March 1997 (revised July 1997) 



Vector Meson Masses in Chiral Perturbation Theory 



J. Bijnens", P. Gosdzinsky'' and P. Talavera** 



"Department of Theoretical Physics, University of Lund 
Solvegatan 14A, S-22362 Lund, Sweden. 

^NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark. 



PACS: 12.39.Fe, 12.40.Yx, 12.15.Ff, 14.40.-n 



Abstract 

We discuss the vector meson masses within the context of Chiral Perturbation Theory 
performing an expansion in terms of the momenta, quark masses and l/N^. We extend 
the previous analysis to include isospin breaking effects and also include up to order p^. 
We discuss vector meson chiral perturbation theory in some detail and present a derivation 
from a relativistic lagrangian. The unknown coefficients are estimated in various ways. We 
also discuss the relevance of electromagnetic corrections and the implications of the present 
calculation for the determination of quark masses. 



1 Introduction 



The relevance of vector meson masses and mixings in the determination of the quark masses is weU 
known. A review is in Virtually all the information used in that paper was at the lowest order 
in Chiral Perturbation Theory and at most the leading nonanalytic corrections are included. The 
main parts of the information used for quark mass ratios were pseudoscalar, baryon and vector 
meson masses and mixings. The full one-loop corrections to the pseudoscalar meson masses and 
mixings was performed in ||] and a recent review on their relevance to quark mass ratios is in . 
The subject of baryon masses in chiral perturbation theory is a subject of a lot of current research. 
For recent results see and references therein. This paper is devoted to a similar study for the 
vector meson masses. 

Chiral Perturbation Theory for Vector Mesons using a heavy meson formalism was first intro- 
duced in Q. We will refer to this as the Heavy Meson Effective Theory (HMET). There the first 
correction of order was evaluated. This was later extended to include isospin breaking effects 
in |6). In that reference an estimate of the electromagnetic effects using a short-long-distance 
matching calculation was also performed. The result for the quark masses from — to mixing 
agreed well with standard expectations, but the K*^-K*^ mass difference did not. The size of 
the chiral corrections of order is rather large so a more complete calculation seemed necessary. 
We therefore perform in this paper the full vector meson mass and mixing corrections to order p'^. 
We also discuss the electromagnetic corrections. 

The paper is organized as follows. In Sect. ^ we give the main notation and conventions. 
In the next section we explain the HMET and list all of the terms necessary for the calculation 
of the vector meson masses. Then in Sect. ^ we explain the connection between the relativistic 
descriptions and the HMET. In Sect. || we then give the model estimates for the various parameters 
in the HMET, followed in Sect. ^ by a description of the calculation of the main result in this 
paper, the vector meson masses to order p*. We then discuss how we include the electromagnetic 
correction in Sect. |^ and the numerical results in Sect. |^. We then give the main conclusions and 
the implications for ratios of quark masses in Sect. I 

In App. ^ we quote the integrals we use in all relevant kinematical domains and in App. ^ 
we give the approximate analytical formulas for the vector masses, as described in Sect. ||. 

2 Basic Ingredients 

In this section we put some of the basic notation. The notation is essentially the same as in 
and§. 

The pseudo-Goldstone boson fields can be written as a 3 (8) 3 special unitary matrix 



U — exp 



iV2U 
F 



(1) 



where 



H 



7r+ K+ 

V6 



K 



(2) 



and F ^ F^^ 92.4 MeV. 

The relativistic vector meson fields are introduced as a 3 (8) 3 nonet matrix (in the Nc limit) 



V2^ V2 ^'^ 



k: 



K*+ 



(3) 



1 



In what follows (A) denotes trace of A. We also use the quantities 

X+ = u^Xu^ + ux^u 
X = 2Bo{s + tp) 



u,. 



-*'fl — fJ,'- 

^{udf,u^ + df,u) (4) 



In (H), we have ignored the vector and axial external sources. We will also set the pseudoscalar 
source, p, equal to zero and s will be the diagonal matrix {mu,md,ms). The parameter Bq is 
given by 

{0\qq\0) ^-F^Boil + Oimg)). (5) 
The transformations under a chiral symmetry transformation x gj^ G SU{3)l x SU{3)n 

are: 



u - 


^ gnUgl 




X ~ 


gnxgl 




u - 


^ gRuh^ = 


hug\ 


R - 


hRh) 


for 


- 


hV^h) + 





(6) 



Eq. (0) also serves as the definition of h, which depends on g^, gn and 11. 



3 Vector Meson Chiral Perturbation Theory 

The Chiral Corrections we are interested in are long distance corrections to a propagating "bare" 
vector mesonQ. We can therefore calculate these corrections in a systematic fashion. We cannot in 
a similar way calculate the chiral corrections to vector meson decays of the kind: V PP since 
that is inherently shorter distance. "Hard" loop corrections to the vector meson masses will be 
included in the constants in the effective lagrangians we will write. An example of this already at 
tree level will be given in Sect. ^. 

It will also be nice to have explicit power counting present in the Lagrangian. We therefore use 
a formalism similar to the Heavy Quark Effective Theory Q and the one used for Baryons||l^ pT| . 
The main part of the mass, the vector meson mass in the chiral limit, we will remove explicitly by 
introducing a velocity v with v'^ — 1 and the chiral mass my. Vector momenta are then treated 
as residual momenta compared to myv. I.e. ky = mvv + p and we only refer to p in the space 
time dependence. We introduce the 3x3 matrix Wf^ for the effective vectors: 



W„ 



V2 



V2 



k: 



K* 



V2 



(7) 



Under the chiral symmetry group VF^ transforms as hW^h^ with h defined in Eq. (||). 

The fields in are to be understood as only containing annihilation operators. The creation 
operators are contained in . In order to have a proper spin one field we impose the condition 



vW = 0. 



(8) 



similar formalism is of course possible for axial vector mesons. The relevance for quark masses is much 
smaller, since masses and mixings are much less well known in this sector 



2 



This condition is enforced by putting in the Lagrangian always the combination P'^"Wa and its 
Hermitian conjugate. Here 

In the remainder this projector should be understood. For the calculation presented in this paper 
it can be forgotten, since the v ■ W component never contributes in any of the diagrams. 
Let us now proceed to construct the relevant Lagrangians. The terms at order are: 

C^^^v{WlW^)+xo{Wl){W''). (10) 

The first term, though allowed by the symmetry, can be removed from the action through the 
choice of the reference momentum. The second term is small because of Zweig's rule, we will 
therefore treat it as the same order as the quark masses, 0{p^). The term at order p is: 

Ci = -i{Wl{v■D)W^')-ixl{Wl){{v■D)W^') + 

ig{{WlWM^pe>^'"'^ + W [{Wl){W,u^) + v^e'^'^^P . (11) 

For the purpose of counting, for constants in the lagrangian that are suppressed by l/A^c, we 
will count them as order p^. So the coefficient xq needs to be kept in the loop diagrams but not 
xi and g' . Similarly, for the p^ Lagrangian below, the 1/Nc suppressed coefficients appearing in 
loops we neglect. 

For the higher orders we can use the lowest order equation of motion to remove terms. So 
we will never encounter {v ■ D)W and {v ■ D)W^ . We will also not explicitly mention terms that 
vanish as the external vector and axial vector fields vanish. These never contribute to the masses. 
In addition, since we treat xi as a small quantity the effects of this term can be reabsorbed in the 
other parameters by the field redefinition 

W'^ = W^ + ^{W^). (12) 

The field still transforms under the chiral group in the same way as W^. In fact the xf 
remainder can also be removed by slightly changing the coefficient in ([l^ ) . 

At order p^ there are several terms. First there are those whose coefficients are fixed by Lorentz- 
invariance. At the level of the HMET Lagrangian these constraints can be directly derived using 
reparametrization invariance[p^. The Lagrangian at order p^ is then 

£2 = ai({Wl, W''}x+) + a2{D^W!D^Wn + asiD^WlD^W) 

+ ai{[{DpWl,W,} ~ {DpW,,Wl}] Uo.)e^'"'^ + a5{{Wl,W''Wu'')+ae{Wlu,W^u'') 
+ a7{{wlW, + Wf^Wl) uf^u") + as{{WlW, + W^W^) u^u^") 

+ ag{Wlu^'W,v}' + Wlu'W.u'') + aio(W^^(v • u)W''iv ■ u)) + aii({W^^, W^} (v ■ uf){13) 
In addition to the terms leading in 1 /Nc of Eq. ( p^ ) we also have 

4 = ai2{{Wl){x+Wn + {x+Wl){Wn)+au{WlW^){x+) 

+ai4Wl){W^){x+)+ai,{D^W!){D^Wn+aie{D^w!){D''W,.) (14) 



contributing to the vector masses directly. The terms of Eqs. (|l0|), ( pT|) and (^3|) are all those 
that will appear within loops in the approximations used here. 
The terms fixed by reparametrization invariance are: 

— 1 g 
a-2 = 7i 04 = 7^ ■ (15) 

We have not performed a full classification of the p^ and terms, this is work in progress, 
but only of those that can contribute to the vector meson masses. We have shown the subleading 
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ones in l/Nc- These are ah the terms needed to absorb the infinities occurring in the present 
calculation. 

There are no terms at tree level of order p'^ that contribute to the vector masses. The argument 
is as follows: In order to have p'^ , we need three derivatives or one insertion of quark masses and 
one derivative. To have an invariant term contributing to the masses we also need a W^, and 
Lorenz invariance requires an extra v. If the v is contracted with the vectors it vanishes because 
V -W ^ Q. If u is contracted with a derivative, it has to act on one of the vectors. This term then 
vanishes by the equations of motion as mentioned above. 

At we have four terms contributing to the masses: 

A = c,{Wlx+W^X+)+C2{{Wl,W^}x+X+) 

+cs{x+{D^W^,,D^W}) + c^ix+iD^WlD^W}) (16) 

The last two terms in fact do not contribute to the masses at order p^. The reason is that the 
external momentum is proportional to a mass difference that is itself of order p^ . So to the masses 
these terms only contribute to order p^. For a similar reason there are no contributions from terms 
with four derivatives. 

The l/Nc suppressed terms are: 

A = C5{x+){x+{WlW^'})+ce{x+Wl){x+Wn 

+cr{x+X+){WlW'^)+cs {{x+X+Wl){Wn + {x+x+W^){Wl)) 
+c,{Wl) {W^) (x+X+) + cio(x+) (x+> {Wlwn 

(X+) {{X+WD {Wn + ix+W^) {Wl))+ ci2 (x+) (X+) {Wl){Wn (17) 



4 Relation to Relativistic Vector Meson Lagrangians 

The usual parametrizations of vector mesons in chiral Lagrangians is done in a relativistic for- 
malism. There are various popular versions of this. A review of some of them exists in p^ . A 
discussion relevant to the relation between the different ways of parametrizing them can be found 
in 01 . The relation in the context of the functional integral approach can be found in ||lj] . 

Let us first treat the case of a single noninteracting Vector Meson. The relevant Lagrangian is: 

C^ = -\v'^-'V^. + ^m'V'^V^, (18) 

with Vp^i, = d^Vu — di,Vfj_. This Lagrangian produces three propagating modes and one constraint 
equation. The latter makes the non relativistic limit a little more tricky than just a naive identi- 
fication. We define first a parallel and perpendicular component of the vector field with respect 
to the velocity v: 

= P^.K + v^{v ■ V) = Vl^ + Vf,V\\ . (19) 

This is similar to the Heavy Quark Effective Theory where a projector {l + v-j)/2 is introduced[||. 
The Lagrangian of ( p^ ) then becomes 



4-- ■ 2 

I- 

2 



+ \{{vd)V\\f-\{d,V\\f . (20) 

We can now split both Vj_ and V\\ in its creation and annihilation parts (keeping in mind the 
projector acting on the effective field): 
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We can then identify my = m + • • • where the dots are a small quantity. We also assume that 
the residual momentum dependence described by W^, VFy and Hermitian conjugates, is small 
compared to mvv. This is where we restrict to the one vector meson sector. 

The only terms in the action that are still proportional to my is (myVFij • W^||)/2. So, in this 
limit we precisely have • = or the constraint we assumed in the previous section v ■ W — Q. 
The symbol W used in the previous section always included the projector P of (^; i.e. it is VFl. 
One could also alternatively integrate out the V\\ component in a functional integral language. 
Here we will just remove it using the equations of motion: 

W\\ = — 9uM^r + \ {dJv ■ d)W'; + W - iv ■ df] W\\ ] . (22) 
my m,y J II J 

This equation can simply be solved iteratively. If applied to the kinetic terms this leads to the 



constraints given in Eq. (15) 



In practice we can choose slightly different values of my to do the reduction to the HMET. 
This is precisely the freedom that leads to the presence of the term proportional to Ay in Eq. 

An alternative procedure, that is in practice somewhat easier to implement, is to use tree level 
matching. We calculate certain processes with the relevant kinematics both in the relativistic 
formulation and in the HMET formulation. We then determine the HMET constants in terms of 
the constants in the relativistic Lagrangian by requiring equality between the two formulations. 
Here one has to watch out somewhat. The equivalence is at the level of iS-matrix elements so 
when comparing Green functions we can only easily compare for the perpendicular components 
and one should remember the extra factors of l/(2my) because the relativistic and HMET fields 
are differently normalized. 



5 Coefficients in the Lagrangian 

The total number of parameters that contribute to order p'^ to the vector meson masses is very 
large. It is obviously too large to be fitted simply from the data on the masses alone. In our 
numerical results we will therefore use several estimates of the parameters. 

The most important ones are the extra 1/Nc assumptions used. Here we add a suppression 
factor of for each factor of l/Nc- The loops themselves are in principle also suppressed by 
factors of l/Nc- But large logarithms and combinatorial factors can in practice make up for this 
extra loop suppression. We therefore still take the loops into account. For the parameters this 
means :0 

g' = ai4 = = ai6 = C5 = • • • = Ci2 = . (23) 

Here equal to zero means that they do not contribute to the masses when the l/Nc expression 
is taken into account, xq should be included inside loops and at tree level, xi is not physically 
relevant as described earlier. ai2 and ai3 only contribute at tree level. ai3 is like a vector meson 
sigma term. On the masses it only contributes just like Ay and is thus of no relevance here. In 
the end there are two possibly relevant l/N^, parameters, Xq and ai2- 

First the vector — vector — pseudoscalar terms. Here we only need g. It was estimated in ||^ 
using the chiral quark model with a value oi g — 0.375. Assuming .g^quark = 1, this changes to 
g — 0.5. In Q the value of g was fixed using a VMD argument for lu ttj, this yielded g = 0.32 
in reasonable agreement with the previous estimate. A double VMD estimate from tt" — > 77 leads 
to essentially the same result. The ENJL prediction for this vertexjl^ is lower g — 0.27. Fitting 
to r(p° — + T"""), r(/0^ ^ 77r+) and T{p ^ tttt), together with various assumptions about VMD 
and the KSRF relation, leads also to values roughly within the above range. We will discuss the 
implications for the vector meson masses of these and other values in Sec. 

^As mentioned earlier, derivatives on external vector field count as for the calculation of the masses. 
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The terms with more derivatives and vector fields only are as described earlier fully determined 
by Lorenz invariance, a2 and a^. 03 never contributes to the masses. 

The vector — vector — pseudoscalar — pseudoscalar interaction terms are more difficult to 
obtain. In Q a good description of Vector meson phenomenology is given by the model III: 

/^iii = + IrnUlv, ~ |r,)'> . (24) 

Here transforms nonlinearly under the chiral group as — s- hV^h^ + hdf^h^ and 

V^. = d^V, - duVf, - ig [F^, K] . (25) 

In this model there are three diagrams possibly contributing to ttV^ scattering. These are depicted 
in Fig. 0. The effects of these have to be described by the pointlike interaction in the HMET. We 
have chosen this model because most properties of Vectors and Axial vectors are well described by 
the Yang-Mills type models. The type "hard" contributions of Fig. ^ cannot be simply described 
using equations of motion as was done in the previous section. For these we use the matching 
procedure. The internal propagator can both be a pseudoscalar or an axial vector. Similarly in 
Fig. |I]c the vector in the vertical line has a small momentum, not a small deviation from myv. 
Diagram Fig. |^ reproduces the pion interactions from the kinetic term and the 02 term of the 
HMET in addition to: ^ 

03 = -a? = as = 75 • (26) 

ZlTiv 

This doesn't contribute to the masses. The diagram of Fig. |^a is reproduced in the HMET via 

= = (27) 



Higher momentum dependences of the nV scattering can be produced by the axial vector 
intermediate states. The V^A^t: vertex in this model can be described by the two terms 

iAi{(v,, - |r^j [u,,D^A,])+iA2{A^ [u,,V,,,]). (28) 

in addition to the kinetic terms for the axial nonet. The field transforms here as R in Eq. (||). 
In the HMET they can be described via 

aio = -2aii = -^^^^^^^-(^1+^2)'. (29) 

The coefficients A\ and A2 should then be chosen to reproduce the rather uncertain a\ width. 

The above model did not contain much internal vector meson properties. As an example of a 
model that does have some internal structure we will use the ENJL model. For a review see [nsl. 



The coefficients needed here were obtained in and 1 15 . In order to obtain the a-i from there we 
have to perform the HMET reduction. We also have to include the diagrams of Fig. |l|, keeping 
in mind that in ^ a different vector representation was used where the vector-pseudoscalar- 
pseudoscalar vertex contains three derivatives. What it in the end corresponds to is that we have 
the above contributions in terms oi g ~ l/{2y/2gv) and Ai = A^^^; A2 — A^^'^ and in addition 
terms coming from the pointlike relativistic contribution of Fig. ^b. These extra terms are 



E 


^(1) 


^(2) 

E _ 
°s " 2mv 


5 — 


4my ' 




^(4) 


x(5) 
E _ (>V 

°9 " 2my 


E 


8 — 


4my ' 



.(3) 
E _ "V 



(30) 
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(a) (b) 



(c) 



Figure 1: The three possible relativistic diagrams contributing to pseudoscalar- vector scattering. 

The symbols gv-, A'^'^\ yl^^^ and ^(i' - '^) are defined and their values in terms of the ENJL param- 
eters given in |T5[ . For the contributions to the masses, ay and ag only appear in the combination 
07 + ag. To obtain the numerical values we can use the ENJL relations: 



2/y 

^(3) ^ ^(2) + + 2.10) « -0.075 . (31) 

JAJV 

The numerical values correspond to fv = 0.20, /a = 0.097 H and Lg + Lio = 1.6 • lO^^Jiel. 

The remaining constants to be estimated are oi, ci and C2. For these we will use an ENJL 
model estimate. The Yang-Mills like models normally just assume these terms to be zero. These 
are the vector meson masses in the large Nc limit. We could in principle determine those from 
the measured meson masses but we will here also produce an ENJL calculation of them. In the 
ENJL model the two-point functions were well described by a VMD picture, see We will 

therefore do the same thing. We will fit the two-point function Iiy^'^\—q^) to Zv/{My — q^). 
We do this fit in the Euclidean region to remove the artifacts due to non-confinement in the ENJL 
model. We observe that a good fit can be obtained with Zy independent of the quark masses and 
will use the fitted My to obtain the coefficients ai, ci and 02- Notice that we fit the squared of 
the mass and not the mass itself, like one should expect in an effective theory. With this we take 
into account the rcsummation of the terms 1/mv suppressed inside these coefficients. The values 
chosen for the ENJL parameters are Gs — 1.216, Gy — 1.263 and A^^, — 1.16 GcV which gives a 
good fit to a large number of hadronic parameters. We have also used the way suggested in Q 
to determine the vector masses. This is another recipe to remove the pole aX = Q in n^^"* . The 
differences are within the quoted errors. We have also done the fit using various ranges of the 
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N 



I 
I 

\ 



Figure 2: The loop diagrams contributing to the masses in the HMET formulation. The dot is 
an L\ vertex. The crosses are £2 vertices. 

quark masses and various ranges of . The results are all within the quoted values. We will use 
the values 

ai = (0.15 ± 0.01) GeV"^ ci = (-0.08 ± 0.03) GeV"^ C2 = (-0.06 ± 0.02) GeV"^ (32) 

The vector mass in the chiral limit obtained from the ENJL model with the above parameters 
is about 0.630 GeV. The model uncertainty on these estimates increases the error on a\ to a 
somewhat larger value. 



6 Determination of the Masses 

The method we will use to determine the masses is to compute the inverse propagator and then 
we look for its zeros. As mentioned earlier we perform the calculation fully in the HMET. We 
are free to choose the external momentum in a way that simplifies the calculation. Our choice 
is that the external momentum is always proportional to the chosen velocity v. This together 
with the constraints v - W = ^ removes a lot of irrelevant contributions. A disadvantage of the 
HMET are that there are more possible diagrams, corresponding to the extra operators with 
fixed coefficients. In the relativistic formulation we only have two possible diagrams. Then we 
have no simple powercounting because the large mass appears in the nonanalytic parts. In the 
nonrelativistic case we have more diagrams. These are depicted in in Fig. ||. 

As argued before we only have to include the effects of the terms up to inside the loops. In 
practice this means that we have to diagonalize the masses including the terms proportional to 
and fli inside the sunrise type diagrams, depicted in the first three diagrams of Fig. |[ We also 
diagonalize the lowest order pseudoscalar mass term to take ir^ — rj mixing into account. 

We proceed in the following way. We define a set of external vector fields with = ^extfj.'^ext 
that corresponds to the usual isospin basis. We then define a 2nd set of "internal" vector fields 
— W'j°(^Tj'^( that diagonalize the lowest order mass terms, i.e. the Ay, ai, xq and Si (see 
next Section) terms. Similarly we define a set of pseudoscalar meson fields tt*^ that diagonalize 
their lowest order mass term, H = tt^^T*^. 
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In terms of these the sunrise diagrams depicted the first three diagrams in Fig. ^ as a function 
of the incoming flavour a and outgoing vector h with residual momentum = PextVfj, can be 
described via: 



M=1.8 c=l,9 



(33) 



The function K is defined as follows: 



m 



M 



V ■ q+ Pext - 



1 



Pext 

mv 



1 



2mv V ■ q + Pext - ^fTla 

(34) 

K can be easily constructed from the integrals given in App. |A|. Arria is the vector meson mass 
minus the chiral limit value, i.e. AttIq = ma — my- to_a/ is the mass of the pseudoscalar. We 
have disregarded in these expressions terms of higher orders found in the integrals, i.e p\^^ and 
Pext^'TTLa, their numerical effects are found to be small and mainly affects de (j). For the we 
actually find a second zero in the inverse propagator for some of the parameter sets. Adding the 
omitted higher order terms in the sunrise diagrams moves this extra zero beyond the range of 
vahdity of CHPT as discussed in Sect. ^. These extra terms correspond to replacing m^^ in the 
last term of K, Eq. (|3|), by m\j + 2pext^ma 



Pext- 



The tadpole contributions, last diagram in Fig. |^ can be given by 



E 



2F2 



f,Mmx„T!:l,}{T'\{T"'\x}}) 



2 

F2 



E 

M=i.a 



- y 



M=l,8 



as 



1 f 



+ as + ail; 



^(2a9 
a 



+ aio) 



(35) 



The traces can be easily performed but lead to rather lengthy expression. In App. ^ we have 
given the expressions in the limit where internal lines are the nonmixed states. The numerical 
difference, especially the part due to the 7r° — 77 mixing, is rather small, less than half an MeV for 
all masses. 

To equations (^) and (^5|) we still have to add all tree level contributions to obtain the 
inverse propagator to order p^. The corrections to the masses of the p+, i4r*°, K*^ are given 
by values of p^xt for which the corresponding amputated two point function vanishes. For the 
/9°, w, the situation is slightly more complicated, because these mix with each other. Here the 
corrections to the masses are given by the values of p^xt for which the determinant of the matrix 
of the amputated two point function vanishes. Once we have found the three solutions, we have to 
identify the particle to which they correspond. This can be done by analyzing the null eigenvectors 
at the value of p^xt where the determinant vanishes. These null eigenvectors define the physical 
basis for the p, oj and 0, and they do not necessarily form an orthogonal set since the physical 
states appear at different values of Pext- 



7 Electromagnetic contributions to the Masses 

The subject of the electromagnetic contributions to masses is in fact very old, an early attempt is 
in In 1^ a matching calculation was performed. We will use the results of that paper here. 
The main electromagnetic effects can be described by the three terms 

h{QWl){QWn+52{[Q,Wl] [Q,W^^])+5^{Q^{WlW^'}), (36) 

with as numerical values Q 

5i = 0.27 ±0.11 MeV ^2 = -1.4 ± 0.4 MeV and cJg = 0.27 ± 0.11 MeV . (37) 
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The numerical results from pC( | cannot be easily compared. The model used there makes rather 
drastic assumptions about the high energy behaviour of the various formfactors needed. The bad 
matching found in Q should however be kept in mind. There are several effects included in pjj , |20| ] 
that in the language of Q are higher order. The main effect of these is the 717 intermediate state. 
The effect of this we model in the present work by including the Si term in the diagonalization 
used in the previous section and by using the measured vector meson masses inside the loops. 

Another possibly large effect is the electromagnetic part of the pseudoscalar meson masses in 
the loops. This is in fact the leading quark mass effect in the electromagnetic corrections Here 
it is included by using the physical pseudoscalar masses within the loopsQ. 



8 Numerical Results 

As a general inputs we choose m„ = 4.273 MeV, nid — 7.727 MeV and iris = 146.0 MeV as well 
as Bq = nr^o/i'iTT'u + iT^d)- These masses are compatible with the small deviation from the Gell- 
Mann-Oakes-Renner relation and the quark mass m — rriu + md determined in l22|. The other 



values are found via the ratios quoted in Q. We also take the coefficients given in Eqs. (31), ( |37| ) 
as fixed parameters. We also fix the ai mass appearing in Eqs. (|2^) via = \/2 mv and choose 
as a subtraction point ^ = 0.77 GeV. That leaves 7 parameters to vary, mv, xq, ai, 012, ci, C2 
and g. We evaluate the masses in different scenarios: 

Scenario /, we use the ENJL estimates for oi, ci, C2, and the set of parameters given in Eq. 
(^). We also use = 6.5 MeV and g = 0.32 with my = mp. We found a large uj(j) mixing 
with a large p^xt dependence. There is also a very large negative correction to the (/)-mass resulting 
in the (f) being the lightest state. Most of these problems are caused by the very large contribution 
from the sunrise diagrams. We can use xo and ai2 to try to solve these problems, would lead 
to a very large ujp^ mass difference so we will use ai2 to cancel the (jjcj) mixing. Requiring the w(j) 
mixing to be small in a reasonable range of Pext leads to 012 ^ 0.16 GeV^^. 

A very similar result can be obtained with g = 0.27 as suggested by the ENJL model. It has 
the masses somewhat better than scenario I but there is no qualitative improvement, we therefore 
do not discuss it further. 

This result is not in good agreement with the experimental data as shown in Table We 
should keep in mind that the values for oi, ci, C2 and 012 we have used are rough estimates, and 
that xq = 6.5 MeV is a tree level value which can receive corrections mainly from the sunrise 
diagram at 0{p^). 

In the remainder we therefore use a least squares fit to determine viable sets of parameters. 
We have minimized the following function: 



mp+—mpo\^ / (mp+ + mpQ)/2\ ^ / (mjf.+ + m^.o)/2 



1 MeV y V 10 MeV / V 20 MeV 

\ 2 / \ 2 / \ 2 



1 MeV / V 1 MeV / V 40 MeV 



(38) 



The masses should be understood here as the difference of the CHPT calculation and the 
measured values. All data are from Ref. |^ except for the — p° mass difference[^. Here 
we have chosen the errors so that the quantities mainly determined by TOj are allowed a larger 
variation. We will then have to check afterwards if the found sets have a reasonable convergence 
in the various orders. 

We will present five different fits, all of which are acceptable. The input parameters can be 
found in Table |^ and the resulting deviation from the observed masses in Table |^. In that table we 
also present the predictions for the various mixings. These are evaluated as the two-point function 
in the isospin basis at the relevant value of the momentum. Notice that in most fits these were 
not input so they present genuine predictions. The ujcj) mixing has a large p^xt dependence. We 

'^The effects on the mixings of pseudoscalars is higher order and we neglect it. 
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Scenario 


I 


II 


III 


IV 


V 


VI 


my (GcV) 


0.7685 


0.84304 


0.85801 


0.84250 


0.74516 


0.75935 


xo (GeV) 


0.0065 


0.00578 


0.00756 


0.00581 


0.0195 


0.0198 


ai (GeV-^) 


0.15 


0.16371 


0.12484 


0.16926 


0.12166 


0.11990 


ci (GeV-^) 


-0.08 


0.13582 


0.15531 


0.13718 


0.17327 


0.17307 


C2 (GeV-^) 


-0.06 


-0.12964 


-0.08406 


-0.13457 


-0.0373 


-0.03873 


ai2 (GeV-^) 


0.16 


0.10284 


0.10070 


0.10279 


-0.01937 


-0.00937 


5 


0.32 


0.25749 


0.26507 


0.25863 


0.11926 


0.14352 



Table 1: The input parameters for the various scenarios. The units are given in the first column 



Scenario 


I 


II 


III 


IV 


V 


VI 


p+ (GeV) 


-0.1666 


-.00001 


0.0001 








K*+ (GeV) 


-0.1658 


+0.0009 


-0.0023 


0.0004 


-0.0003 




K*'> (GeV) 


-0.1708 


+0.0009 


-0.0025 


0.0004 


-0.0003 




(GeV) 


-0.1665 


-0.0005 


0.0001 


-0.0005 


-0.0001 




uj (GeV) 


-0.1638 




0.0001 


-0.0005 


-0.0001 




(GeV) 


-0.3177 




0.0003 


-0.0009 


+0.0003 


+0.0001 


wp"(p") (GeV) 


-0.0023 


-0.0024 


-0.0016 


-0.0025 


-0.0016 


-0.0014 


p> (0) (GeV) 


+0.0004 


+0.0009 


+0.0008 


+0.0008 


+0.0006 


+0.0005 


uj(j) (u;) (GeV) 


-0.035 


-0.042 


-0.048 


-0.043 


-0.071 


-0.072 


oj(f> {(p) (GeV) 


-0.056 


+0.015 


+0.009 


+0.011 


-0.055 


-0.050 




18 


0.49 


0.19 


0.52 


0.06 


0.017 



Table 2: The deviation from the observed vector masses, for the various scenarios. Also given 
is the total standard deviation, this V^mm defined in Eq. (^8|). The units are given in the first 
column. A dash means agreement to better than 4 digits. 
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have given it both at the w-mass and the 0-niass. The pcf) mixing is more stable and we quote it 
only at the 0-mass. For definiteness we have quoted the p^uj mixing at the p'^-mass. Using the 
w-mass instead is not visible with the precision quoted. 

Scenario II is a minimum found with a fairly large value of g. All masses can be fit within 
the errors assumed in Eq. (|38|). All parameters are of the expected order of magnitude here, the 
main change is the flip in sign of ci. As mentioned earlier, there is a second spurious zero in 
the determinant of the inverse two-point functions of the neutral states about 60 MeV above the 
(^-mass. In order to check the sensitivity to higher orders of this phenomenon and the possible 
variations of the parameters we have also performed a fit with the higher order parts of the sunrise 
diagrams included. This pushes the spurious zero to about 230 MeV above the 0-mass and varies 
the input parameters within less than one standard deviation from the previous fit. That case is 
in the tables as scenario III. 

A much better minimum can be obtained for a value of g around 0.1. This is quoted in scenario 
V and VI. The latter is again with the higher order parts of the sunrise diagrams included. 

As can be easily seen we obtain reasonable agreement with the observed values for all the 
mixings as derived in These were: 

p°LO = -2.5 MeV; = 0.3 MeV and uj(j){(l)) = (8 to 14) MeV . (39) 

We should keep in mind that the latter two have possibly large final state corrections. 

The agreement with the p^uj mixing in scenario III is not so good. We have therefore added 
to Eq. dH) also 

/ p°^(p°) + 0.00254 GeV 
V 0.0002 GeV 

Including this then gave a minimum within the errors expected of the input parameters. This we 
have given as scenario IV. 

We have not discussed the p, dependence in the two-point function because this can be absorbed 
in a change in the parameters in principle if we include the 1 /A^c-suppressed terms as well. In 
practice we have to fit the parameters anyway. We have checked that putting p = 0.9 GeV instead 
we can get as good as a fit as scenario III and IV with similar changes in the input parameters. 

To see the convergence of the series for the various quantities we show in Table |3| the various 
contributions corresponding to scenario IV. They are ordered by the order of the contributions. 
The sunrise diagrams contain pieces of order and higher. The 02 term is a term with two 
derivatives but as explained earlier these count as for the calculation of the vector masses. For 
the masses of the neutral states we have quoted the various contributions to the unmixed inverse 
two-point function at the correct mass. The effect of the mixing on the masses can be seen by 
comparing the total line to the real masses as obtained in Table ^ for scenario IV. As can be seen 
the convergence for most quantities is acceptable but the (/)-mass is very slowly converging. 




9 Conclusions 

We have in this paper extended the calculation of the vector meson masses in CHPT beyond p^g 
and included also the isospin breaking effects. The problem that the corrections break several 
observed regularities can be repaired by the introduction of the higher orders. The large breaking 
of the equal space relation, — mx* = rriK* — rUp, and the large lu — 4> mixing are also solved. 
In addition the problem of the K*^ — K*^ mass difference^ with usual quark mass ratios can be 
resolved. The parameters we obtain from a fit to the masses are very good with reasonable choices 
of all the parameters. The orders of magnitude of the various parameters expressed in GeV are all 
natural with the possible exception of xq in some scenarios. In addition the predictions we have 
for all the mixings using these fitted parameters are also very good. 

We conclude that the vector meson masses and mixings can be well described within GHPT 
with the standard quark mass ratios. 

^In Ref. |20| the latter problem was discussed at order including the electromagnetic corrections and a 
different solution proposed. 



12 







K*+ 






LiJ 


4> 






UJ(p{(j)) 


ai 


12.3 


154.5 


158.0 


12.3 


12.3 


300.1 


-3.55 








xo 














11.6 


5.8 








8.2 


Sun 


-140.8 


-126.2 


-127.2 


-140.8 


-143.6 


-329.5 


3.26 


2.6 


-90.3 


ai-tad 


26.3 


53.2 


53.6 


26.3 


26.3 


80.5 


-0.44 








Oj-tad 


30.1 


70.7 


71.7 


30.1 


29.0 


58.1 


-0.28 


0.4 


-41.1 


ai2 














15.0 


182.3 


-2.16 


-1.5 


134.2 


Cl 


0.2 


3.2 


5.7 


0.2 


0.2 


107.8 


-0.10 








C2 


-0.4 


-105.9 


-106.1 


-0.4 


-0.4 


-211.6 


0.20 








02 


-3.2 


-1.5 


-1.7 


-3.3 


-2.2 


-18.4 











em 


1.5 


1.5 


0.06 


1.5 


0.3 


0.4 


0.54 


-0.6 


-0.2 


Total 


-74.0 


49.5 


54.0 


-74.0 


-51.4 


175.5 


-2.53 


0.8 


10.8 



Table 3: The various contributions to masses and mixings for scenario IV. All units here are 
MeV. The contributions are labeled by the coefficient in front of the vertex for the tree level 
contributions. The loop diagrams are the sunrise diagrams (Sun), the tadpoles with the ai vertex 
(ai)-tad and the other tadpoles (oi-tad). The electromagnetic corrections from Q (em) and the 
sum are in the final lines. 



Note added 

After this paper was submitted we became aware of a paper where the value of g has been 
determined from the rates of r ^ lotti^. The values obtained there, are in good qualitative 
agreement with those used here, they correspond to g « 0.29 and g w 0.32 depending on the 
number of experimental bins included. 
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A Integrals 



We quote here the integrals we use. First the simple scalar ones 

d'^q l(g2) _ im2(m'*) r 



A-log(^ 



(27r)2 g2 _ ^2 iQ^2 

with A = l/e — 7 + log(47r) + 1 and d = 4 — 2e. The others can all be expressed in terms of iJo 

d'^q 1 1 



(41) 



iJo = 



(27r)2 V ■ q — uj + irj q^ — + irj 



167r2 

The function F{uj, m) is 



'2luX + uj 



F{ijj,m) = 



2V^ 



:l0g 



^2 + 21og(^ 



UJ + \/ up- ~ 



w — vtiJ — mr 



4(m^ - u?)F{uj, m) 



(lo > m) 



(42) 
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— Lo^ 



77 UJ . 

o - arctg^=== (-m < w < m) 



' log , , , + 2j7r < -m) . (43) 



We now define 

(27r)2 V ■ q — Lv + ir] q^ — w? + 177 

i [JiW^; J2g^v + JiV^V,,; JiiQu^Va + ff^a?;^ + Si^aW;:.) + hv^V^V^, JsQ^v + JlV^V^] .(44) 

These can all be expressed in terms of Jq and /x^: 

(45) 

Jo - f^Pm) (46) 

(47) 

(48) 

(49) 
(50) 
(51) 

The expressions for those with higher powers of the "heavy" propagator can be derived by taking 
derivatives of the above expressions with respect to oj. The only difficulty there is the derivative 
of the function F. This is given by 

^^ = ^(-i^(-,-)-l). (52) 

B Approximate Analytical Expressions 

In this appendix we collect the explicit formulae for the masses, in the approximation of non- 
diagonal fields inside the sunrise diagrams. 
For the tree level contribution we find: 

= 4aiBo(m„ + rud) + a2Pext + ^QBlmurridCi + 16Bo(m^ + m^)c2 - ^2 + g^s 

5p° = 4alBo(m„+md)+a2pLt + 8^o("^« + "^d)(cl+2c2) + y + -^3 

5 S 

5p^oj = 4aiBo(m„-md)+8S^(m^-m^)(ci+2c2) + -^ + ^ 

D 6 



Jl 


= ojJq + 






1 

d- 1 


((m^- 


h 


= LoJi — 


J2 


Ji 


= CJ J2 + 


vn? 


h 


= W J3 - 


2J4 


Je 


= J2 




J7 


= TO^Js 





,0, 



Si 

3V2 

i-Bo 

_6i_ 
9V2 



Su = 4aiBo{mu + md) + a2pixt + ^^^{ml + ml){ci + 2c2) + % + ^63 
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S 2 

6(j) = SaiBorUs + a2pl^t + 16Bo™?(ci + 2c2) + + 

5 

5K*~^ = 4aiBo{mu + ms) + a2plxt + '^QB^mumgCi + 16BQ{ml + rnl)c2 - 62 + -S3 

2 

6K*° = AaiBo{md + m^) + a2Pext + 16-Bo"^dmsCi + 16Bo(m^ + ml)c2 + -S3 
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For the tad-pole type diagrams, we find: 



5uj(j) 



5p^ = aiBo— ^iJ.^o{muC^n + mdc^) + 2/z^+(to„ + md) + iiK+{mu + rris) + HK^imd + ms) 

+ ^(^2e.,-(.2 + ^)e,-(c^ + ^)e.o 

+ iiniruusl + mds^) j - ( + S^fo + 2E^+ + (c^ + y )S^o + (s^ + _)l]^ 



.(^2e,+ + (s2 + ^)e, + (c2 + ^)e,o 



2b2 



If 4 

bK*^ = aiBo^iiJ,^o{mucl + -mss'^) + ii„+{mu + md) + 2ij,K+{mu + ms) + iJ-Ko{md + rns) 
+ Unimus; + -msC^)J - (^4S^+ + 2I]ko + 2S„+ + (— + — + s^)E^ 



,5s2 2.SC 



+ - = + c^)E,.) + ((^ + ,c)e, + (ii - ,c)6„) 



2 / 49 

6K*° = ai Bo -p;2 ( Mtt" ("^dCp + - ms s ) + (m„ + m^) + (m„ + ms ) + 2/ij^o (m^ + TOs ) 

9 4 9A bi / ^ „ _ ,5c^ 2sc 

+ Hr,{mdsl, + -msc^) j - — ( 2Sk+ + 4E^o + 2S,+ + (_-_ + s^)E^ 

469 / \ 

2 5c)e^ + (i_+sc)e„,j (53) 



We have introduced the following notation: 
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07 + 08 + an ae ag oio 

b,=a, + 62 = y + - + — (54) 



and 

,2 



16^ 7^ SP - mp (^/ip + ^ 805 + 2(07 + 08 + on) 



o 2 / , 209 + Oio \ , ■ 

6p = nip A^p + -I „ 2 Q ] 7: (55) 



MP = -TTJT^log-^ Ep=mp(^/iP 

mp 2o9 + Oio 
167r2 806 + 4o9 + 2oio J 

We have also defined : 

s s c c . . 

c — cos 9, s = sinfc', Cp = c~\ — ^, Cm = c ■=,Sp = s-\ — Sm = s ^, (56) 

v3 v3 v3 v3 

For the sunrise type diagrams, we find: 

Sp+ = A|^c2[,y,p+] + [fc+,A-*°] + [fcO,if*°]+2[7r+,H + ^s2[^o,p+] 
Sp° = A |2 .2 [rj, u;] + lc' [ry, p"] + [fc+, + K*'] + 2 [tt", + [tt", p^] | 

,5pV = V2A{[A:+,A-*+] - [fc",X*0]} 

= A |^c2[r7,c.] + 2s'[rj,p°] + [fc+,if*+] + [fc", X*"] + 4^+, P+] + ls'[n°,u;] + 2c2[^o,pO] 
= V2A{[k+,K*+] + [k",K*'']} 
6^ = Ai^^c^[r,,^] + 2[k+,K*+] + 2[k^,K*"] + ^s'[7T", 

r 1 rs 1 1 



SK*'> = 



A {y[7?,X*°] + ^hK*'] + + [k\p+] + + [fc°,, 

We have introduced 

A.|!, p.Ql.(l + ^-||^)j,(™„,™„-™„-.^..) (57) 

Where the [G, fi] function is related with the one defined in (]3^). The physical 77^ and ttq are 
related to the isospin eigenfields 778 and ttq by 

rjg — cos9 rjp — sin0 ttq ttq = sinS rjp + cos0 ttq, (58) 
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